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Abstract Monte-Carlo generator with pho- 
ton jets radiation in collinear regions for the pro- 
cess e + e~ — > 77 is described in detail. Radiative 
corrections in the first order of a are treated ex- 
actly. Large leading logarithmic corrections com- 
ing from collinear regions are taken into account 
in all orders of a by applying the Structure Func- 
tion approach. Theoretical precision of the cross 
section with radiative corrections is estimated to 
be 0.2%. This process is considered as an ad- 
ditional tool to measure luminosity in forthcom- 
ing experiments with the CMD-3 detector at the 
e+e- collider VEPP-2000. 

1. Introduction 

For the first time the process 

e+(p+) + e"0_) -> j(ki) + 7(fc 2 ), 

was considered in the classical papers by L.Brown 
and R.Feynman [1], I.Harris and L.Brown [2] in 
early 1950s and then revised in 1973 by F.Berends 
and R.Gastmans [3]. Due to the large magni- 
tude of the cross section of this process, it can 
be exploited as an independent way to measure 
luminosity. Precise determination of luminosity 
is a key ingredient in all experiments which study 
hadronic cross sections at e + e~ colliders. As a 
rule, a systematic error of luminosity measure- 
ments represents one of the largest sources of un- 
certainty which can cause significant reduction of 
an accuracy of the hadronic cross sections normal- 
ized to luminosity. The process of two-gamma- 
quantum annihilation has essential advantages for 
luminosity measurements with respect to those 
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based on events of Bhabha scattering or annihi- 
lation into a muon pair. Indeed, the cross section 
value estimated for large angles is of the same 
order as that of Bhabha scattering. Events of 
this process have two collinear photons at large 
angles providing a clean signature for their selec- 
tion among other detected particles. The CLEO- 
II collaboration was the first to show in practice 
how the combined application of the processes 
e + e~ — > e+er , and 77 helped to achieve a 

1% accuracy of the luminosity measurement [J]. 

It is worth noting that the dependence of the 
Born cross section on the photon polar angle 9 is 
not so steep as that of Bhabha scattering and it 
is symmetric under the transformation 6 — > n — 
facilitating a study of the systematic uncertain- 
ties on the detector acceptance. In addition, it is 
free of difficulties related to both radiation and 
Coulomb interaction of the final state particles. 
It is also of utmost importance that correspond- 
ing Feynman graphs do not contain photon prop- 
agators affected by vacuum polarization effects, 
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Figure 1. Feynman diagrams for two-gamma- 
quantum annihilation. 
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Therefore, knowledge of this cross section with 
radiative corrections (RC) at the level of per mill 
accuracy is urgently needed. On the other hand, 
it is a purely quantum electrodynamics (QED) 
process giving large background while studying 
hadronic processes with neutral particles in the 
final state. These reasons are the main motiva- 
tion to consider this process with precise radiative 
corrections and create a Monte-Carlo generator to 
simulate e + e~ — > 77 events. 

2. Cross section of the process e + e — » 77 
with the first order a corrections 

The differential Born cross section for the two- 
gamma-quantum annihilation in the Born ap- 
proximation reads 



that they are not detected and its value does not 
exceed some small quantity Ae <C £■ 

Consider now the process of three-gamma- 
quantum annihilation which can be treated as a 
radiative correction for two-gamma-quantum an- 
nihilation: 

e+( P+ ) + e~(p_) -> 7(k x ) + 7(fca) + 7(^3). 

For the first time the analytic expression for this 
process was obtained by M.V. Terentjev in [6]. 
A much simpler way to obtain the same expres- 
sion was suggested in [7] , based on the chiral am- 
plitude method, when all three hard photons are 
emitted in the so called non-collinear region (out- 
side narrow cones). The cross section is given by 
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where s = (p + +P-) 2 = 4e 2 , e is the beam energy 
and (3 — v/c = y/l — ml/e 2 , c\ = cos6>i, 9\ = 
k\p—. It is assumed that both final photons are 
detected and their polar angles with respect to 
the beam directions are not too small (6*12 ^> 
m e /e). 

Following the well known results obtained in [5] 
and considering the RC due to emission of virtual 
and soft real photons, this part of the cross sec- 
tion can be written as 



da 



B+S+V 

+K SV 

3 1 + Cl 
2 1 - ci 

•ln 2 i^l 



dan I 1 



(L - 1) 2 In 



Ksv = 



3 2(1 



In 



I -ex 



1 



Ae 

e 

1 + 
1 1 



1 - ci 

1 + ci + 2 1 - ci 



+ (ci -> -Ci] 



(2) 



where L = ln(s/m 2 ) is a large logarithm. For 
s ~ 1 GeV 2 , L 1=3 15. The energy of the radiated 
soft photons is implied to be sufficiently low so 
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where Xi = hp-, Xj_^ hp+, i=l,2,3, x { = k°/e, 
at = cos(6i), 9i =p_ki. 

Energy-momentum conservation allows to de- 
termine kinematics of the final photons: 



Xl + x 2 



X3 

Xl 
C2 : 



X1C1 + X 2 C 2 + X3C3 
1 — X3 







1 • 2 ib ■ 

1 — X3 sin 
cos 2 1 + (1 



\2 ■ 2 ib 

x 3 ysm f 



• 2 ib 

X3 sm 



X1C1 + X3C3 
X2 



i> = fcife 3 . (4) 



The sum of this cross section with those describ- 
ing soft and virtual photons radiation does not de- 
pend on inner parameters. It allows to construct 
a MC event generator to simulate three-photon 
events and to take into account proper selection 
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criteria of a given experiment as well as include 
specific detector imperfections. 

2.1. Matching NLO and higher-order (HO) 
corrections 

It is known that photon jets are radiated in 
collinear regions along the motion of electrons and 
positrons give the dominant contribution to the 
cross section. So, in order to achieve a theoreti- 
cal precision of about per mill, all enhanced HO 
corrections must be taken into account and com- 
bined with NLO corrections. The opening angle 
of these narrow cones is small and obeys the re- 
strictions: I/7 <C 6q -C 1. As a rule, its value is 
chosen as 9q — 1/^/7. Since the photon radiation 
outside these cones is not enhanced, it is sufficient 
to consider only radiation of one photon at large 
angles [9 to keep the theoretical accuracy at the 
per mill level. 

For completeness, the cross section with one 
hard photon emission in the collinear region is 
presented below. It can be obtained using the 
method of quasi-real electrons [8] similarly to [9] . 
The theorem of factorization [10] of hard and soft 
photons permits to treat RC in the leading log- 
arithmic approximation in all orders of pertur- 
bation theory in terms of Structure Functions of 
electron and positron, T>(z,L). This fact allows 
one to write the differential cross section for the 
process of two gamma-quantum annihilation in- 
tegrated inside the collinear region as follows: 
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where z = 1 — x is electron (positron) energy after 
photon emission with energy x. 

The shifted Born cross section with reduced en- 
ergies of the incoming electrons and positrons has 
the form 
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where the scattering angle 9\ is given for the orig- 
inal cm. reference frame of the colliding beams, 
z\ and Z2 are the energy fractions of electron and 
positron just before collision after radiation pho- 
ton jets. When z\ and Z2 tend to unity, this cross 
section is transformed to (fT]). One can see that a 
part of this cross section has a term proportional 
to large logarithm L = ln(s/m 2 ), due to collinear 
photon emission, and it is already contained in 
Structure Functions. Therefore, to match NLO 
and HO corrections and exclude double counting, 
the term proportional to large logarithm must be 
removed. The remaining non-leading correction 
referred to as a compensator should be combined 
with the cross section describing three-photon 
production to cancel the dependence on the aux- 
iliary parameter 9q. 

3. Total cross section of the process e + e — » 

77 +n7 

In the following, we summarize the main fea- 
tures of the matching procedure as implemented 
in the code MCGPJ [9]. Adding the higher-order 
RC in the leading logarithmic approximation to 
the complete one- loop result (NLO), the master 
formula for the resulting cross section can be rep- 
resented as follows: 



da e 



77+( " 7) = / AzxD{zx) / dz 2 V(z 2 ) ■ 



■da (zi,z 2 ) (l H Ksvj Q(cuts) + 



dx 



7T J X 
A 



x" \ . On x 2 ' 
z H In — H 



da (z,l) +da (l,z) 
4a 3 



Q(cuts) 

4(i + 4) 



X;>A 

w-e >8i>9 



^(l-c 2 )(l-c|) 



-two cyclic permutations 



drO(cuts), (7) 



where T>(z) is the smoothed representation for the 
Structure Functions according to [12] . A factor 
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1/3 in the last term takes into account the iden- 
tity of the final photons. The all variables are 
defined above. 

This expression contains the logarithmically 
enhanced contributions due to emission of pho- 
tons at all powers of a in collinear regions and, 
as it will be shown later, provides the cross section 
accuracy of about 10~ 3 . The first term describes 
radiation of photon jets which is approximated by 
the convolution of the Structure Functions with 
the shifted Born cross section (s' = SZ1Z2). The 
step functions Q(cuts) stand for particular cuts 
applied. The sum of the last two terms provides 
cancelation of the auxiliary parameters A and 8 . 
A detailed comparison was performed between 
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Figure 2. Cross section dependence on the auxil- 
iary parameter Ae, y/s = 900 MeV. 



the results obtained with MCGPJ and the MC 
generator based on [TT] for the cuts modeling the 
CMD-2 event selection criteria at the cm. energy 
y/s = 900 MeV. These cuts are: 

• The polar angles of the two most energetic 
photons must be inside a range 1.1 < #i ; 2 < 



1.1 



• Acollincarity must obey \6i + 9 2 — tt\ < 0.25 
and l^i - 02 1 -tt\< 0.15 

• The energies of the two most energetic pho- 
tons must be larger than half of beam en- 
ergy. 
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Figure 3. Cross section dependence on the auxil- 
iary parameter do, y/s — 900 MeV. 



In Fig. [2] the cross section as a function of the 
Ae is shown when other parameters are fixed ac- 
cording to selection criteria pointed above. It can 
be seen that there is a broad plateau, where the 
cross section deviations lie within a band with a 
width of ~ 0.2%, whereas Ae runs by more than 
two orders of magnitude. Only for large values of 
Ae some trend appears which can be explained 
by the omitted terms proportional to Ae/e. Sim- 
ilar dependence is seen in Fig. [3j where the cross 
section variations with the auxiliary parameter 
00 are presented. Only for extremely small values 
of 9q, when a condition 1/7 <C #0 is not valid, 
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the cross section falls down by ~ 0.1% only. It 
is an important result, which confirms that the 
cross section does not depend on the auxiliary 
parameters Ae and 9q. It is important to re- 
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Figure 5. Relative cross section difference calcu- 
lated with the MCGPJ code and generator based 
on Ref. [5]. 



Figure 4. Relative cross section difference calcu- 
lated with the MCGPJ code and generator based 
on Ref. [5] versus the angle |A0|. 



liably estimate the total theoretical precision of 
this cross section with RC. In order to quantify 
the theoretical accuracy, an independent compar- 
ison has been performed with the MC event gen- 
erator based on , where only first-order a cor- 
rections are treated exactly. It was found that the 
relative difference of the cross sections is larger 
than 1% for small angles A9 < 0.1 radians and it 
is about ~ 0.6% for an acollinearity angle ~ 0.25 
radians. The simulation results are presented in 
Fig. |4j It is seen that the difference practically 
does not depend on the choice of the value Ad 
with accuracy ±0.1%. The same difference was 
studied as a function of beam energy when the 
acollinearity angle A8 was fixed to 0.25 radians. 
The results are shown in Fig. [SJ In this case the 



difference slowly increases with energy from 0.3% 
to 0.7%. It is an important fact, which means 
that for this energy band the radiation of two and 
more photons (jets) in the collinear region con- 
tributes to the cross section by amount ~ 0.5% 
only. As the uncertainty of this correction is at 
least a few times smaller than the magnitude of 
the correction itself, we can conclude: a theoret- 
ical precision of the cross section with RC is cer- 
tainly better than ~ 0.2%. 

4. Conclusion 

The main features of the Monte-Carlo genera- 
tor to simulate events of the process e + e~ — > 77 
are described. The theoretical precision of the 
cross section with RC is ~ 0.2%. It was achieved 
due to the application of the Structure Function 
approach which allows one to match the enhanced 
contributions (HO), coming from the collinear re- 
gions, with the cross section containing NLO cor- 
rections. It is proposed to use the generator as a 
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complementary tool to measure the collider lumi- 
nosity compared to the regularly used processes 
of Bhabha scattering or e + e~ — > fi + \i~ . The pre- 
cision of collider luminosity determination repre- 
sents one of the largest sources of systematic er- 
rors which can cause significant reduction of the 
accuracy of hadronic cross sections, which are, as 
a rule, normalized to luminosity. The considered 
cross section is rather big, events have a clean 
signature in the detector and can be easily rec- 
ognized and selected. From the theoretical point 
of view, the main advantage of using the process 
e + e" — > 77 compared to others is the following: 
the cross section does not contain the first order 
of a corrections due to the vacuum polarization 
effects and there is no FSR and Coulomb interac- 
tion of particles. 
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